Defect Solitons in Parity-Time Symmetric Optical Lattices with Nonlocal Nonlinearity 
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The existence and stability of defect solitons in parity-time (PT) symmetric optical lattices with 
nonlocal nonlinearity are reported. It is found that nonlocality can expand the stability region of 
defect solitons. For positive or zero defects, fundamental and dipole solitons can exist stably in the 
semi-infinite gap and the first gap, respectively. For negative defects, fundamental solitons can be 
stable in both the semi-infinite gap and the first gap, whereas dipole solitons are unstable in the first 
gap. There exist a maximum degree of nonlocal nonlinearity, above which the fundamental solitons 
in the semi-infinite gap and the dipole solitons in the first gap do not exist for negative defects. 
The influence of the imaginary part of the PT-symmetric potentials on soliton stability is given. 
When the modulation depth of the PT-symmetric lattices is small, defect solitons can be stable for 
positive and zero defects, even if the PT-symmetric potential is above the phase transition point. 

PACS numbers: 42.25.Bs, 42.65.Tg, 11.30.Er 
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I. INTRODUCTION 

The physical systems with parity-time (PT) symme- 
try have attracted much more attention in recent years 
[lHlOJ. Bender et al. found that PT-symmetric Hamil- 
tonians can have entirely real spectra although these 
Hamiltonians are non-Hermitian |l|, |2J|. In optics, PT- 
symmetric structures can be constructed by inclusion of 
gain or loss regions into waveguides, which make the 
complex refractive-index distribution obeying the condi- 
tion n(x) — n*(~x) |3j-|5(. Many unusual features stem- 
ming from the PT-symmetry have been found, such as 
power oscillation [f| |l0|, LLL| , absorption enhanced trans- 
mission '9], nonlinear switching structure [12|, [I3], and 
unidirectional invisibility [14J . In experiments, the PT- 
symmetry breaking in complex optical potentials was ob- 
served firstly by Guo et al in a double- well structure fab- 
ricated through a multilayer Al x Gai^ x As heterostruc- 
ture with varying concentrations [9]. Riiter et al observed 
the non-reciprocal wave propagation in an active PT- 
symmetric coupled waveguide system based on Fe-doped 
LiNbO^, in which optical gain is provided through two- 
wave mixing using the photorefractive nonlinearity [10j. 
Those theoretical and experimental results led to the pro- 
posal of a new class of PT-symmetric synthetic materials 
with intriguing and unexpected properties that rely on 
the non-reciprocal light propagation [lij, (Hf ■ 

In optics, nonlinearities in the PT-symmetric systems 
have been considered by many authors [3, 13; ToL espe- 
cially in the PT-symmetric optical lattices 4], y, [U | , and 
some new kinds of soliton were found and investigated 
0) fl7H20| . Defect solitons can be formed when a local 
defect is introduced into the optical lattices, which have 
been widely studied in the regimes of photonic crystals, 
waveguide arrays, and optically induced photonic lattices 
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[21H23J . Due to its unique properties based on the defect 
guiding phenomena which can be controlled through vari- 
ation of the defect parameters [24| , defect solitons have 
the potential applications for the all-optical switches and 
routing of optical signals [25| . Zhou et al have studied the 
defect modes in the PT-symmetric optical lattices (18J . 
Defect solitons in PT-symmetric optical lattices and su- 
perlattices with local nonlinearity have been studied, and 
stable solitons are found mainly in the semi-infinite gap 
01. 



It is noteworthy that the nonlinearity in the photore- 
fractive media, in which Riiter et al. have observed the 
non-reciprocal wave propagation [10J, is nonlocal due to 
the diffusion mechanism of charge carriers [26- 28] . The 
nonlocality of nonlinear response exists in many real 
physical systems, such as photorefractive crystals 129- 
[28j , nematic liquid crystals [29l - [31j , lead glasses [32|, [33| , 
etc. The nonlocality can drastically modify the proper- 
ties of solitons and improve the stability of solitons [3J- 
[361 ]. Therefore it is worthy to study the properties of 
solitons in PT-symmetric optical lattices with nonlocal 
nonlinearity 



In this paper, we study defect solitons in PT- 
symmetric optical lattices with nonlocal nonlinearity 
The properties of nonlocal defect solitons are quite dif- 
ferent from those in local media, and the nonlocality can 
expand the stability regions of soliton, especially in the 
first gap. It is found that the stability of nonlocal de- 
fect solitons depends on the defect, the degree of nonlo- 
cality, and the PT-symmetric potentials. The influence 
of the imaginary part of the PT-symmetric potentials is 
studied, and the example of stable soliton in the PT- 
symmetric potentials above the phase transition point is 
given. 



II. THEORETICAL MODEL 

We consider the propagation of light beam in the PT- 
symmetric defective lattices with Kerr-type nonlocal non- 
linearity. The evolution of complex amplitude U of the 
light fields can be described by following dimensionless 
nonlinear Schrodinger equation (NLSE), 



.dU ld 2 U 



dz 2 dx 2 



+ p[V(x) + iW(x)}U + nU = 0, (1) 



where x and z are the transverse and longitudinal co- 
ordinates, respectively, and p is the depth of the PT- 
symmetric potentials. V{x) and W(x) are the real and 
imaginary parts of the PT-symmetric potentials, respec- 
tively, which are assumed in this paper as 



V(x) = cos 2 (x)[l + eexp(-a; 8 /128)], 
W(x) = W sm(2x). 



(2) 
(3) 



Here e represents the strength of the defect, which is ex- 
pressed as a super- Gaussian profile [l{|. The parameter 
Wo represents the strength of the imaginary part of the 
PT-symmetric potentials compared with the real part. 
The nonlinear refractive-index change n satisfies 



44 = ^i 2 > 

ox z 



(4) 



where d stands for the nonlocality degree of the nonlinear 
response. This type of nonlinear response with a finite re- 
gion of nonlocality exists in many real physical systems, 
for instance, all diffusion- type nonlinearity [28l I37I. l38j . 
orientational-type nonlinearity [3l|, [39j , and the general 
quadratic nonlinearity describing parametric interaction 
|40l |41| . Although the diffusion nonlinearity in photore- 
fractive media is anisotropic, Eq. (j4|) can still be used 
to describe its nonlocal feature theoretically when we ig- 
nore the anisotropic effect [431. Moreover when d —¥ 0, 
the above equation reduces to n — \U\ 2 , and Eq. ([T|) re- 
duces to that for local case. When d ^ 0, Eqs. ([T]) and 
(HJ denote a nonlocal NLSE with an exponential-decay 
type nonlocal response. 

The PT-symmetric lattices governed by Eq. ^ have 
Bloch band structures when e = 0. The band diagram 
can be entirely real when the system is operated below 
the phase transition point (W§ = 0.5) [3, Q- The Bloch 
band structure obtained by the plane wave expansion 
method for p = 4 is shown in FigQJa). From Fig. QTa), 
we can see that the region of the first gap decreases with 
increasing Wo- The first gap disappears when Wo = 0.5, 
which is the phase transition point. When the system is 
above the phase transition point, i.e. Wo > 0.5, the band 
structure becomes complex. Figures (Hb)-[TJd) show the 
profiles of the PT-symmetric potentials for p = 4 and 
Wo = 0.1 with positive defects (e — 0.5), zero defects 
(e = 0), and negative defects (e = —0.5), respectively. 

We search for stationary solutions to Eqs. (JIJ and (@| 
in the form U — f(x)exp(ibz), where f(x) is a complex 
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FIG. 1. (color online) (a) Band structure for the PT- 
symmetric lattice with p = 4; The blue regions are bands 
while the blank regions are gaps, (b)-(d) Lattice profiles for 
lattices with p = 4, Wo = 0.1, and (b) e — 0.5, (c) e = 0, and 
(d) e = —0.5, respectively. Solid blue and dotted red lines 
represent the real and imaginary parts, respectively. 



function satisfies equations, 



bf 



2dx 2 



-p[V(x)+iW(x)]f + nf, 

»-45-i/i". 



(5) 
(6) 



The solutions of defect solitons are gotten numerically 
by the modified squared-operator method [43| from Eqs. 
§5§ and © and shown in the next section. To elucidate 
the stability of defect solitons, we search for perturbed 
solution to Eqs. ([T]) and (|4]) in the form U(x, z) = [f(x) + 
u(x,z) + iv(x,z)]exp(ibz), where the real [u(a;, z)] and 
imaginary [v(x, z)] parts of the perturbation can grow 
with a complex rate 5 upon propagation. Linearization 
of Eq. ([T]) around the stationary solution f(x) yields the 
eigenvalue problem, 

1 d 2 u 

Sv = - -— - + (n — b)u + p(Vu - Wv) 

2 ax z 

+ Re[f(x)} [°° 2G(x-0u(0Re[f(0]dt 

J — oo 

+ Re[f(x)] f°° 2G(x - £)/m[/(£)M£)d£, (7) 

J —oo 



Su 



ld 2 v 



(n — b)v — p(Wu + Vv) 



2dx 2 
Im{f(x)} / 26'!, -i)M/K>HiKl; 



— OO 



-Im[f(x)] 2G(x-ORe[f(OHOdt (8) 

J — OO 

Here G(x) = [l/(2d 1 / 2 )] exp(-|a;|/d 1 / 2 ) is the response 
function of nonlocal nonlinearity. Above eigenvalue prob- 



lem is solved numerically to find the maximum value of 
Re{5). If Re(6) > 0, solitons are unstable. Otherwise, 
they are stable. 



III. NONLOCAL DEFECT SOLITONS 

In the PT-symmetric defective lattices with nonlocal 
nonlinearity, we find two types of defect solitons for pos- 
itive, zero, and negative defects, respectively. The first 
type is the nodeless fundamental solitons, which can ex- 
ist stably in the semi-infinite gap for all kinds of defects 
or in the first gap for negative defects. The other type of 
defect solitons, which exist in the first gap for all kinds 
of defects, is called dipole solitons in this paper, because 
most of them have two significant intensity peaks. 

For positive defects, we assume e = 0.5 and the results 
are shown in Figs. [5] - 2) Figures [5Ja) and [U[b) show 
that the power of soliton [defined as P = f_°° \f(x)\ 2 dx] 
for both fundamental and dipole solitons increases al- 
most linearly as increasing of propagation constant. As 
the nonlocality degree d increasing, the power of soliton 
increases too. The positive defect solitons vanish when 
the propagation constant is below the cutoff point, whose 
value does not depend on the nonlocality degree. This 
feature is similar to the case of traditional uniform lat- 
tices in nonlocal media 1351] . 




FIG. 2. (color online) Soliton power versus propagation con- 
stant for positive defect with (a) Wo = 0.1 and (b) Wo — 0.15. 
(c) The unstable growth rate Re(5) for soliton solutions in (b). 
For all cases e = 0.5 and p = 4. 

When Wq = 0.1, the linear stability analysis shows 
that both fundamental and dipole solitons are stable in 
their whole regimes of existence in the semi-infinite gap 
and the first gap, respectively. For comparison, dipole 
solitons in the PT-symmetric defective lattices with local 
nonlinearity are stable only in a small region near the 
edge of the first gap |1||. It shows that the nonlocality 
expands the stability region of soliton. As Wo increases, 



the stability region of defect soliton decreases. Figure 
&[c) shows the maximum growth rate Re(8) for solitons 
with Wo = 0.15. We can see that both fundamental 
and dipole solitons are unstable for the large propagation 
constant, and their stability regions become very narrow. 
We find that all fundamental solitons are unstable and 
their stability region disappears when Wq > 0.24. For 
the dipole solitons in the first gap, the maximum value 
for existence of stable soliton is Wq — 0.16. 
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FIG. 3. (color online) The complex fields (solid blue: real 
part; dotted red: imaginary part) and the refractive-index 
changes (dashed green) for soliton solutions at (a) 6 = 6, 
d = 1; (b) b = 6, d = 3; (c) b = 2.5, d = 1; (d) b = 2.5, 
d = 3; (e) 6 = 6, d = 1; and (f) b — 2.5, d = 3; respectively. 
For all cases e = 0.5, p = 4. For (a)-(d) W = 0.1 and (e)-(f) 
Wo = 0.15. 

The field profiles of the positive defect solitons are 
shown in Figs. [3][a)-[3Jf) for different Wo and nonlocal de- 
grees d, which correspond to the circle symbols in Figs. 
[5£a) and HJb) . Figures OUa) and OUb) show fundamen- 
tal solitons in the semi- infinite gap for Wq =0.1, while 
dipole solitons in the first gap are shown in Figs. Efc) and 
[3jd). One can see that the poles of dipole solitons are 
located inside the central channel of the lattice. As the 
nonlocality degree increases, the amplitudes of solitons 
increase but their shapes change very little. Figures [3Je) 
and[3^f) show the defect solitons for Wo = 0.15. As the 
Wo increases, the imaginary part of the field increases. 
The propagations corresponding to those solitons in Figs. 
E£a)-[3Jf) are shown in Figs. IUa)-|3Jf). The evolutions are 
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FIG. 6. (color online) The complex fields (solid blue: real 
part; dotted red: imaginary part) and the refractive-index 
changes (dashed green) for solitons in the semi-infinite gap at 
(a) b = 4, d = 1; (b) b — 4, d — 3; or solitons in the first gap 
at (c) 6 = 1, d — 1; and (d) 6 = 1, d = 3, respectively. For all 
cases e = 0, p — 4 and Wo = 0.1. 



FIG. 4. (color online) (a)-(f) Evolutions of defect solitons 
corresponding to those in Figs. |3ja)-[3jf), respectively. 



simulated based on Eqs. ((T|) and (j4]), and 1% random- 
noise perturbations are added into the initial input to 
verify the results of linear stability analysis. We can see 
that both fundamental solitons and dipole solitons are 
stable for Wq = 0.1 but unstable for Wq = 0.15. 
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FIG. 5. (color online) Soliton power versus propagation con- 
stant for zero defect with (a) Wo = 0.05 and (b) Wo = 0.1, 
respectively, (c) The unstable growth rate Re(5) for soliton 
solutions in (b). For all cases e = and p = 4. 



Figures [5] and |5] show the results for zero defects 
(e = 0). The properties of defect solitons with zero de- 
fects are similar to those with positive defects. From 
Figs. [H^a) and[5]Jb), we can see that the cutoff points 
approach the lower edges of gaps. For Wo = 0.05, funda- 
mental solitons are stable in the whole semi-infinite gap, 
whereas dipole solitons are stable in the whole first gap. 
When Wo =0.1, dipole solitons are unstable for the large 
propagation constant, and their stability region becomes 
narrow, as shown in Fig. EJc). The stability region of 
dipole solitons vanishes when Wq > 0.14. Fundamen- 
tal solitons are still stable in the whole semi-infinite gap 
for Wo = 0.1, and their stability region begins to reduce 
when W > 0.12. When W > 0.23, all fundamental 
solitons become unstable. Figure [5] shows the complex 
profiles of solitons for zero defects, which correspond to 
the circle symbols in Fig. [5jb). Although the examples 
of solitons in the first gap [Figs. [BJc) andlHJd)] have more 
intensity peaks, the most of solitons in the first gap still 
have two significant peaks. The propagations of those 
examples for Wq — 0.1 in Fig. [6] are all stable as the 
expectation by the linear stability analysis. 

Next we study defect solitons in the nonlocal PT- 
symmetric lattices with negative defects, and results are 
shown in Figs. [7J- [TUJ Fundamental solitons are found 
both in the semi-infinite gap [Fig. [7] and the first gap 
[Fig. [8], whereas dipole solitons exist in the first gap 
[Fig. [7j. It is interesting that there exists a maximum 
nonlocality degree d max for some negative defect solitons. 
When the degree of nonlocal nonlinearity is above the 
maximum value, the fundamental solitons in the semi- 
infinite gap and the dipole solitons in the first gap do no 
longer exist. The maximum nonlocality degree depends 
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FIG. 7. (color online) (a) The maximum nonlocality de- 
gree d max versus the depth of negative defect for fundamen- 
tal solitons in the semi-infinite gap (dashed line) and dipole 
solitons in the first gap (solid red line), (b) and (c) Soliton 
power versus propagation constant for negative defects with 
(b) Wo = 0.05 and (c) W = 0.1, respectively, (d) The un- 
stable growth rate Re(S) for soliton solutions in (c). For all 
cases p — 4 and for (b)-(d) e = —0.5. 




FIG. 8. (color online) Soliton power versus propagation 
constant for fundamental solitons in the first gap with (a) 
Wo = 0.1, e = -0.5; (b) W = 0.1, e = -0.2; and (c) 
Wo = 0.2, e = —0.5, respectively, (d) The unstable growth 
rate Re(8) for soliton solutions in (c). For all case p — 4. 



on the depth of negative defects, as shown in Fig. EJa). 
We can see that d max decreases with increasing of the 
defect depth |e|. We also find similar phenomenon on the 
traditional defect solitons in non-PT-symmetric lattices 
with negative defects. 

For negative defects, the linear stability analysis shows 
that fundamental solitons in the semi-infinite gap are sta- 
ble in their whole existence region for Wq = 0.05. As Wq 
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FIG. 9. (color online) The complex fields (solid blue: 
real part; dotted red: imaginary part) and refractive-index 
changes (dashed green) for soliton solutions in the semi- 
infinite gap at (a) b = 2.8, d = 0.5; (b) b = 2.8, d = 1; 
or soliton solutions in the first gap at (c) b = 2, d — 0.5; (d) 
b = 2, d = 1; (e) b = 2, d = 0.5; and (f) 6 = 2, d = 1, 
respectively. For all cases e = —0.5, p — 4 and Wo = 0.1. 



increases, the stability region reduces. For Wq = 0.1, 
Fig. [T^c) shows that the fundamental solitons in the 
semi-infinite gap are unstable in the most existence re- 
gion except at the edge of the gap when d — 1 , whereas 
they are all unstable when d = 0.5. For comparison, the 
local defect solitons in the semi-infinite gap are unstable 
near the edge of the gap [19( . All fundamental solitons in 
the semi- infinite gap are unstable when Wq > 0.13. The 
dipole solitons in the first gap are always unstable in the 
whole existence region for any value of Wq. 

Differing from the positive and zero defects, fundamen- 
tal solitons can also exist stably in the first gap. Figures 
[Sfa) -IHJc) show the power of fundamental solitons versus 
propagation constant in the first gap. We can see that 
there exists a cutoff point of propagation constant above 
which the fundamental solitons can exist. This feature is 
similar to the case of the positive and zero defects. The 
cutoff point shifts toward the low propagation constant 
with increasing of the depth of negative defect and the 
value of Wo , and the existence region of fundamental soli- 
tons in the first gap increases too. When Wq — 0.1, the 
fundamental solitons in the first gap are stable in their 
whole existence region, and their stability region begins 
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FIG. 10. (color online) (a)-(f) Evolutions of defect solitons 
corresponding to those in Figs. |9ja)-[9jf) respectively. 



to reduce when W > 0.16. When W = 0.2, funda- 
mental solitons is unstable in the most existence region 
except at the edge of the gap, as shown in Fig. [HJd). 
Their stability region vanishes when Wo > 0.23. 

The examples for the negative defect solitons are shown 
in Fig. [9j including the fundamental solitons in the semi- 
infinite gap [Figs. I^a) and[SJb)], the dipole solitons in 
the first gap [Figs. [9jc) andHJd)], and the fundamental 
solitons in the first gap [Figs. H^e) andJHJf)], which corre- 
spond to the circle symbols in Fig. [7£c) and Fig. HJa), re- 
spectively. Their corresponding propagations are shown 
in Figs. HDtal-HOff). Figures EJa) and ITOTa) show the 
unstable fundamental soliton nearby the edge of the gap 
for d = 0.5, whereas , the fundamental soliton nearby the 
edge of the gap for d = 1 is stable as shown in Figs. |9jb) 
and ITUT b) . One can see that dipole solitons in the first 
gap are unstable [Figs. [TUlx) and HUTd)], whereas fun- 
damental solitons in the first gap are stable [Figs. HUT e) 
andfTOTf)] for W a = 0.1. 

In our study, all defect solitons are unstable for large 
Wq when p = 4. However, the solitons in Ref. [4| is 
stable when p = 1 and Wq = 0.45, which is close to the 
phase transition point (W§ — 0.5). So, we finally study 
the solitons with small modulation strength of the PT- 
symmetric potentials. We find that when the modulation 
depth of the PT-symmetric potentials is small, defect 
solitons can be stable for the positive defect and zero 
defect though Wq is above the phase transition point. 
Figures [TlTa)- lriT c) show the examples of defect soli- 
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FIG. 11. (color online) The complex fields (solid blue: real 
part; dotted red: imaginary part) and the refractive-index 
changes (dashed green) of fundamental solitons with (a) pos- 
itive defect e — 0.5, (b) zero defect e = 0, and (c) nega- 
tive defect e = —0.5. (d)-(f) Evolutions of defect solitons 
corresponding to those in (a)-(c), respectively. For all cases 
p = 0.1, Wo = 0.55, b = 0.2 and d = 1. 



tons above the phase transition point, i.e. p = 0.1 and 
Wq = 0.55, with positive defects (e = 0.5), zero defects 
(e = 0), and negative defects (e = —0.5), respectively. 
Figures [TTl d)- ITlT f) show the corresponding propagations 
of those solitons. With the supporting of the nonlocal 
nonlinearity (d — 1), defect solitons with positive and 
zero defects are stable when Wo is above the phase transi- 
tion point, whereas they are unstable for negative defects. 
Here the self-guiding effect by the nonlocal nonlinearity 
plays more important role than the defect guiding effect, 
in balance with the diffraction and the transverse energy 
flowing induced by the imaginary part of Pt-symmetric 
lattice. In an other word, the nonlocal nonlinearity en- 
hanced the real part of the PT potential, thus the nonlin- 
earity can change the phase transition point equivalently 
4]- ' 



IV. SUMMARY 

In conclusion, we have studied the existence and sta- 
bility of defect solitons supported by PT-symmetric non- 
local optical lattices. It is found that the nonlocality 
expands the stability region of defect solitons, especially 



in the first gap. Fundamental solitons can exist stably in 
the semi-infinite gap for all kinds of defects or in the first 
gap for negative defects, whereas dipole solitons exist in 
the first gap for all kinds of defects. Most of defect soli- 
tons are stable for small value of Wo, and their stability 
regions reduce and finally vanish as increasing Wo- As 
an exception, dipole solitons with negative defects are al- 
ways unstable for any value of Wo . It is interesting that 
there exists a maximum nonlocality degree d max for some 
negative defect solitons. We also find defect solitons with 
positive and zero defects can be stable even if the PT- 
symmctric potentials is above the phase transition point 



when the modulation depth of the PT-symmetric lattices 
is small. These properties of the defect solitons in PT- 
symmetric nonlocal lattices are obviously different from 
those in the local PT-symmetric lattices. The various 
types of solitons may provide alternative methods in po- 
tential applications of synthetic PT-symmetric systems. 
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